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(Dated: August 7, 2008) 

Assuming antiferromagnetic orbital correlations to model the pseudogap state in the underdoped 
high-temperature superconductors, we study how this correlation is distinguished from the d-wave 
superconductivity correlation with including the finite-range antiferromagnetic correlation effect. In 
spite of the fact that both correlations have the same d-wave symmetry, the contributions from each 
correlation is clearly distinguished in the spectral weight and the density of states. 

PACS numbers: 71.18.-|-y, 74.20.-z, 74.72.-h 



In the high-temperature superconductors, the origin 
of the pseudogap state is still in controversial. Con- 
trary to conventional superconductors, behaviors asso- 
ciated with opening some gap are observed in various 
experiments above the transition temperature to the su- 
perconducting state, Tc.[l[ Concerning the origin of the 
pseudogap, mainly there are two pictures. One is to as- 
sume that the pseudogap is associated with precursor of 
superconductivity. The other is to assume another gap 
formation which is different from the superconducting 
gap-iiiig 

Recent angle-resolved photoemission spectroscopy 
measurements seem to support two gap scenario. |a] 
In the normal state, the full Fermi surface is clearly ob- 
served. Along this Fermi surface, the gap has been mea- 
sured in the superconducting state. Around the nodal 
region, the gap is well fitted by a d-wave superconduct- 
ing gap. As moving away from the node, the gap deviates 
from the line expected from the simple d-wave supercon- 
ducing gap. The doping dependence of the gap around 
the node and that around the antinodal region are in- 
vestigated. It is found that the doping dependences of 
the gaps in these different regions are qualitatively differ- 
ent. The gap around the node increases as the hole dop- 
ing concentration is increased in the underdoped regime. 
By contrast, the gap around the antinodal region de- 
creases. Furthermore, the gap around the antinode shows 
negligible change across Tc whereas the gap around the 
node shows significant variation. This result suggests 
that those gaps have different origins. This is consis- 
tent with the recent Raman scattering experiments. [lo| 
In the Raman scattering experiments the nodal region 
and the antinodal region are distinguished by symmetry. 
From Big and B2g spectra, the Raman shifts show qual- 
itatively different behaviors. The shift associated with 
the node (B2g) increases as the doping concentration is 
increased while the shift associated with the antinode 
(Big) decreases. 

The two gap picture is also supported by scanning 



tunneling microscopy (STM). About a decade ago, some 
STM measurements supported that the pseudogap and 
the superconducting gap had the same origin because the 
pseudogap seemed to smoothly evolve into the supercon- 
ducting gap as decreasing temperature. However, in the 
recent STM measurements [lli it is found that there are 
two components: One is inhomogeneous and the other is 
homogeneous. The inhomogeneous component is associ- 
ated with relatively high-energy features and it is found 
that that component does not show clear temperature 
dependence around Tc. Boyer et al. extracted the ho- 
mogeneous component from the raw data by taking a 
normalization using the spectra taken above Tc.[ll| The 
gap in the homogeneous component clearly vanishes at 
Tc. By contrast, the inhomogeneous component does not 
show qualitative difference at Tc- It is also found that 
the gap value estimated for the homogeneous component 
is reduced from that estimated from the raw data. It is 
known that the latter gap Ap is scaled by the pseudo- 
gap temperature, T* with 2Ap ~ 4.3fcBT*. This scaling 
relation is consistent with the orbital antiferromagnetic 
correlation with d-wave symmetry, [l^ 

If both of superconductivity and the pseudogap are 
characterized by d-wave symmetry, then the question is 
how we can distinguish these correlations in the physical 
quantities. To answer this question, we study a simple 
phenomenological model for the two gap state for high- 
temperature superconductors. We assume the d-wave 
BCS gap for superconductivity and d-wave orbital an- 
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tiferromagnetic correlation [^, Il2l Il3l . Il4| to model the 
pseudogap. Based on this model we compute the spec- 
tral weight along the underlying Fermi surface and the 
density of states with including the short-range antifer- 
romagnetic correlation effect. In spite of the fact that 
the two gaps have the same symmetry, we find that 
these components are clearly distinguished in the Bril- 
louin zone and in the density of states. 

The model Hamiltonian is given by 
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where Ek - 
2t2 (cos 2kx 



-2tQ (cos fca; + COS k 



COS 2fc, 



At I COS fca; COS fcy — 

Afe — ^ {cos kx — COS ky) , and 

2 (cosfcj; — cosfcj,). The summation with re- 
spect to k is taken over the reduced BriUouin zone, \kx\ < 
TT and \ky \ < tt. The wave vector Q is fixed to Q = (tt, tt) 
for the pure mean field state. In order to include the 
short-range antiferromagnetic correlation effect approxi- 
mately, we assume that Q obeys a Lorentzian probability 
distribution whose peak is at (tt, tt) with broadening fac- 
tor of Caf- We follow the analysis described in Ref.fisj 
in investigating the short-range antiferromagnetic corre- 
lation effect. 

The energy dispersions are obtained by diagonalizing 
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Here 6^+^ = £i±£^ and e^"^ = £i^£i^±« . For the cal- 
culation of the spectral weight, we compute the Mat- 
subara Green's function, Gk] (t) = — (TrCk] {t) cj,| (0)^ 
Transforming to the Matsubara frequency from the imag- 
inary time T, we have 



Gk] (iuJn) = ^ 



Uk1,a («7t) 



(3) 



The Unitary matrix Uk-\,a is computed numerically diag- 
onalizing Eq. IT]). The index a is for the energy bands. 
We have four bands for Ajf^ ^ and Aq ^ 0. 

In Figdl the spectral weight at the Fermi energy is 
shown. The solid line represents the under-lying Fermi 
surface which is given by Ek — fJ-- The chemical potential 
fi is determined from 
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with Hk ^ 2 + <^4-^Cfe|^-{^c_fexcLfe^^ + (^4+QTCfc+QT 

(^C-k-Qic'_k-Qi^ y ^ the doped hole concentration and 
Ns the number of the lattice sites. The expectation val- 
ues are calculated by diagonalizing ^ . We take the hop- 
ping parameters as to = 1, ti = —0.25, and t2 = 0.10 
through out the paper. 




FIG. 1: Spectral weight at the Fermi energy averaged over 
probability distribution with respect to Q. The antiferromag- 
netic correlation length is ^af ~ 4. The solid line represents 
the underlying Fermi surface with Aq^' — and Ao = 0. 



Although we assume a finite antiferromagnetic cor- 
relation length, the mean field calculation shows that 
2Ao/fcBr* ~ 4.8 for the interaction F = 2 and ^af = 4 
on 50 X 50 lattice. Therefore, we may assume a finite 
value of Ao in spite of the lack of the orbital antiferro- 
magnetic long range ordering. 

Figure shows the spectral weight along the Fermi 
surface at fixed Q = (7r,7r). In order to identify the con- 
tributions from the superconductivity correlation and the 
orbital antiferromagnetic correlation, we compare the re- 
sult with or without each correlation. One can see that 
the structure around the superconducting gap node is 
mainly determined by the d-wave superconductivity cor- 
relation. Meanwhile the antinode structure is mainly de- 
termined by the orbital antiferromagnetic correlation. 

Figure [3] shows the dispersion energies along the un- 
derlying Fermi surface. By comparing with Fig[21 we see 
that a part of the energy dispersions is invisible in the 
spectral weight because of the coherence factor effect. As 
shown in Fig. HI the energy dispersion around the node is 
well described by d-wave gap function. Around the anti- 
node, the energy dispersion is also well described by d- 
wave symmetry. However, this component is dominated 
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FIG. 2: Spectral weight along the underlying Fermi surface. 
The vertical axis represents the energy and the horizontal 
axis represents the angle along the underlying Fermi surface: 
(a)Ao/to = 1.0 and Af'^ = 0.3, (b)Ao/to = 0.0 and Af*^ = 
0.3, and (c)AoAo = 1.0 and Af^ = 0.0. The angle is defined 
so that the value of 45 is at the nodal point and the values of 
and 90 are at the anti-nodal points. 
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FIG. 3: The energy dispersion along the underlying Fermi 
surface. The horizontal axis represents angle along the un- 
derlying Fermi surface. 

by the orbital antiferromagnetic correlation as suggested 
from Figl^l This analysis suggests that clear informa- 
tion about the superconducting gap is extracted near the 
node. In fact, from the linear function fitting of the gap 
near the node we find Ag*-^ from FigUl By contrast, 
the gap near the antinode contain both contributions of 
the superconductivity correlation and the orbital antifer- 
romagnetic correlation. It is hard to decompose these 
components even in our idealized model. 

Next we compute the density of states as follows 

For numerical computations, we replace the delta func- 
tion with a Lorentzian. The broadening factor is chosen 
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FIG. 4: The dispersion energy versus | cosfc^, — cosfcy|/2. 




FIG. 5: Energy versus density of states averaged over prob- 
ability distribution with respect to Q for different values of 

as 0.10. In Fig. [5] we show for different values of 

S^AF- We can distinguish each band contribution. The 
peaks around ~ ±1 are associated with the orbital an- 
tiferromagnetic correlation. While the shoulders around 
uj ~ ±0.3 are associated with the superconductivity cor- 
relation. These features become sharp for large i^^i? . But 
the broadening effect due to Q fluctuations are similar as 
shown in FiglS) This result suggests that inhomogeneity 
of the relatively high energy component is not associated 
with the antiferromagnetic correlation effect. 

Figure [7] shows the intersection of D(lS) at = 0.079. 
The same result is obtained at = —0.079. For large 
i^AF, we see some weight outside of the reduced Brillouin 
zone. By contrast, this weight disappears for short ^af- 
As a result, we see a banana shape as seen in STM ex- 
periments. 

To summarize, we study the spectral weight and the 
density of states in the presence of both of d-wave su- 
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FIG. 6: Energy versus density of states averaged over proba- 
bility distribution with respect to Q at ^af ~ 20. Error bars 
are due to probability distribution of Q. 
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FIG. 7: Density of states computed by Eq. ([5]) at a; = 0.079 in 
a quadrant of the Brillouin zone: (a)^AF ~ 50 and (b)f af = 
5. For both panels, we take Aq = 1 and A^^^^ — 0.3. 



perconductivity correlation and the orbital antiferromag- 
netic correlation with taking into account the finite-range 
of the antiferromagnetic correlation effect. Although 
these correlations are assumed to have the same d-wave 
symmetry, it is demonstrated that these components are 



clearly distinguished in the spectral weight and the den- 
sity of states. As for asymmetry of the density of states, 
we need to include impurity scattering effect as demon- 
strated in Ref.lldl. 
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